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Abstract 

The estimates of the gluon concentrations in the classical SU(3)-Yang-Mills field 
modelling confinement are given for the first three states of charmonium whose 
spectrum is tuned by calculating electromagnetic transitions among the mentioned 
levels in dipole approximation. For comparison the corresponding estimates for the 
photon concentration in the ground state of positronium (parapositronium and or- 
thopositronium) are also adduced. 

Key words: quantum chromo dynamics, confinement, charmonium 
PACS: 12.38.-t, 13.20.Gd, 14.40.Gx 



1 Introduction and preliminary remarks 

As was discussed in Refs. [1-3], the natural way of building meson spectroscopy 
might be based on the exact solutions of the SU(3)-Yang-Mills equations 
modelling quark confinement, the so-called confining solutions. Such solutions 
are supposed to contain only the components of the SU(3)-field which are 
Coulomb-like or linear in r, the distance between quarks. In Ref. [1] a num- 
ber of such solutions have been obtained and the corresponding spectrum 
of Dirac equation describing the relativistic bound states in those confining 
SU(3)-Yang-Mills fields has been analysed. It should be noted that the given 
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approach is the direct consequence of the relativistic QCD (quantum chromo- 
dynamics) Lagrangian since the mentioned Yang-Mills and Dirac equations 
are derived just from the latter one. 

Further in Refs. [2,3] the results obtained were successfully applied to the 
description of the quarkonia spectra (charmonium and bottomonium). In its 
turn, the mentioned description points out the linear confinement to be (clas- 
sically) governed by the colour magnetic field linear in r. 

The results of Refs. [1-3] suggest the following mechanism of confinement to 
occur within the framework of QCD (at any rate, for mesons and quarkonia). 
The gluon exchange between quarks is realized in such a way that at large dis- 
tances it leads to the confining SU(3)-field which may be considered classically 
(the gluon concentration becomes huge and gluons form the boson condensate 
- a classical field) and is a nonperturbative solution of the SU(3)-Yang-Mills 
equations. Under the circumstances mesons are the relativistic bound states 
described by the corresponding wave functions - nonperturbative solutions of 
the Dirac equation in this confining SU(3)-field [1-3]. For each meson there 
exists its own set of real constants (for more details see below) aj, Aj,bj, Bj 
parametrizing the mentioned confining gluon field (the gluon condensate) and 
the corresponding wave functions while the latter ones also depend on /xq, the 
reduced mass of the current masses of quarks forming meson [1-3]. It is clear 
that constants a,j, Aj,bj, Bj, /i should be extracted from experimental data 
and such a program has been just realized in Refs. [2,3] for quarkonia. 

The aim of the present paper is to estimate the above gluon concentrations 
for the case of charmonium but we shall restrict ourselves to the first three 
states of charmonium. For comparison we shall also adduce the corresponding 
estimates for the photon concentrations in the ground state of positronium 
(parapositronium and orthopositronium) to emphasize some analogy between 
QCD and QED (quantum electrodynamics). 

Further we shall deal with the metric of the flat Minkowski spacetime M that 
we write down (using the ordinary set of local spherical coordinates r, ip for 
the spatial part) in the form 

ds 2 = g^dxP <g> dx v = dt 2 - dr 2 - r 2 (d$ 2 + sin 2 $d^ 2 ) . (1) 

Besides, we have \5\ = \ det(g^ u )\ = (r 2 sini9) 2 and < r < oo, < < 7r, 
< up < 2tt. 

Throughout the paper we employ the system of units with H = c = 1, unless 
explicitly stated otherwise. Finally, we shall denote L/2(F) the set of the mod- 
ulo square integrable complex functions on any manifold F furnished with an 
integration measure while L%(F) will be the n-fold direct product of L 2 (F) 
endowed with the obvious scalar product. 
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If A = A^dx^ = A c ^\ c dx^ is a SU(3)-connection in the (trivial) three-dimensional 
bundle £ over the Minkowski spacetime, where A c are the known Gell-Mann 
matrices, then we are interested in the confining solutions A of the SU(3)- 
Yang-Mills equations 

d* F = g(*F A A — A A *F) (2) 

with the exterior differential d = d t dt + d r dr + dgdd + d^dip in coordinates 
t, r, (f, while the curvature matrix (field strentgh) for the ^-bundle is F = 
dA+gAAA = F^ v \ a dx^ Adx u and * means the Hodge star operator conforming 
to metric (1), g is a gauge coupling constant. 

For the aims of the given paper we shall use the confining solution of Ref. [1] 
in the form 

Al+^A* t = -^+A U -Al+^Al = ^±^-(A 1+ A 2 ), --L4? = -^+A 2 , 
V3 r * v / 3 r V3 r 

K+^ A l = bir+B^ -A%+±=A% = -(6 1 +6 2 )r-( J B 1+J B 2 ), ~-j=A% = b 2 r+B 2 

(3) 

with all other A c ^ = 0, where real constants aj, Aj,bj, Bj parametrize the 
solution, and we wrote down the solution in the combinations that are just 
needed to insert into the Dirac equation 

V^f = no® (4) 

with the relativistic wave function \P = \& 2 , \& 3 ) of the quarkonium, where 
the four-dimensional spinors \Pj represent the j-th colour component of the 
quarkonium, /x is a mass parameter and one can consider it to be the reduced 
relativistic mass which is equal, e. g. for quarkonia, to half the current mass of 
quarks forming a quarkonium, while the coordinate r stands for the distance 
between quarks. The explicit form of operator T> is not needed here and can 
be found in Refs. [1-3]. 

Finally, for the necessary estimates we shall employ the T 00 -component of the 
energy-momentum tensor for a SU (3)- Yang-Mills field 

T "" = ~L ("^ F ^ ^ + F ^^ 7 V) • (5) 



2 Relativistic spectrum of charmonium 



From the adduced form it is clear that the solution (3) is a configuration de- 
scribing the electric Coulomb-like colour field (components A t ) and the mag- 
netic colour field linear in r (components A v ) so the solution can be rewritten 
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as 

A\ = -(ai + a 2 /2)/r + A 1 + A 2 /2, A t 8 = (a 2 /r - ^2)^/2 , 

A% = (h + b 2 /2)r + B,+ B 2 /2, A% = -{b 2 r + B 2 )V?>/2 . (6) 

As was shown in Ref. [1], after inserting the above confining solution into Eq. 
(4), it admits the solutions of the form (with Pauli matrix a±) 

*, = e*W ),i = l,2,3 (7) 



with the 2D eigenspinor $j = ^ ^ of the euclidean Dirac operator on the 

unit sphere S 2 . The general explicit form of $j is not needed here and can be 
found in Refs. [4]. For the purpose of the present paper we shall adduce the 
necessary spinors below. Spinors % form an orthonormal basis in L 2 (S> 2 ) and 
they can be subject to the normalization condition 

J <Z>}<Z> j( m = J J (l^il 2 + |$ i2 | 2 ) wnWdtp = 1 , (8) 

§2 00 
where (t) stands for hermitian conjugation. 

The energy spectrum e of quarkonium is given by the relation e = lo\ + uj 2 + u 3 
with 



A -A x g 2 a l b l + (n x + a l )J(n\ + 2n l a l + A\)iil + g 2 b 2 1 (n\ + 2n l a^ 
Ul = Ul ( ni , h, A l} = gA 1+ n? + 2n 1 a 1+ Af 

(9) 

uj 2 = uj 2 (n 2 , l 2 , A 2 ) = -g(Ai + A 2 )+ 



-A 2 c/ 2 (a! + a 2 )(6i + 62) - (n 2 + a 2 )V(™2 + 2n 2 a 2 + A|)^g + # 2 (&i + & 2 ) 2 (™2 + 2n 2 a 2 ) 



n 2 . + 2n 2 a 2 + A 2 



2 ~r ^"-2^2 -r ii 2 

fio) 



^3 = ^3(^3, ^3, A 3 ) = gA 2 + 



-A 3 g 2 a 2 b 2 + (n 3 + a 3 )\/( n 3 + 2ra 3 a 3 + A§)//§ + # 2 & 2 (n| + 2n 3 a 3 ) 



n§ + 2n 3 a 3 + A§ 

(11) 

where A x = A x - , A 2 = A 2 + g{B 1 + B 2 ) , A 3 = A 3 - #5 2 , = 0, 1, 2, 
while Aj = ±(lj + 1) are the eigenvalues of euclidean Dirac operator on unit 
sphere with lj = 0, 1, 2, ... Besides 

ai = y/Al - g 2 a\ ,a 2 = ^A 2 - g \ ai + a 2 ) 2 , a 3 = ^/A§ - # 2 a! . (12) 
Further, the radial parts of (7) are given at rij = by 

F n = ([,!>/■ , " (l - ^ , F, 2 = /f^r'S " (l + ^) ? (13) 
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Table 1 

Gauge coupling constant, mass parameter and parameters of the confining SU(3)- 
connection for charmonium. 



g Mo a\ a 2 b\ b 2 Bi B 2 

(GeV) (GeV) (GeV) 



0.46900 0.62500 2.21104 -0.751317 20.2395 -12.6317 6.89659 6.89659 



while at rij > by 



Y 



PjQj 



^-^)L^(r j ) + ^r j Ll a i L + 1 1 (r j ) 



F j2 = iC j Q j r a 'e 




2a, +1, 



J-^Un— 1 



with the Laguerre polynomials Lf^.(rj) 



; = 2f3jr, where, for j - 
Hi - (wi - gA^ 2 + g 2 b\, P 1 = gb x + (3 X , Q 1 = // + wi - gA u Y l 



(14) 
1, Pi = 

[<XlPl ~ 

ga 1 (u 1 -gA 1 )+ga 1 b 1 ]Q 1 +g 2 a 1 b 1 P 1 , Z 1 = [(\ 1 -gB 1 )P 1 +ga 1 /i )]Q 1 +g 2 a 1 b 1 P 1 , 
while for j = 2 one should replace a±, Ai,bi, Bi — > — (a± + a 2 ),—(Ai + 
A 2 ), — (bi+b 2 ), —(B 1 +B 2 ) and for j = 3 one should substitute a±, Ai, b\, B 1 for 
a 2 ,A 2 ,b 2 , B 2 to obtain the corresponding quantities Pj, Qj, Yj, Zj, so that, for 

instance, (3 2 ■ 



I 111 - [uj 2 + g{A l + A 2 )} 2 + g 2 ^ + b 2 ) 2 , f3 3 = y ^ - ^ 3 - g^ 2 
Also it should be noted that the quantum numbers rij of (9)— (11) are defined 
by the relations (for more details see Ref. [1]) 



gA 2 ) 2 + g 2 b 2 . 



ni = 



gb 1 _Z 1 - L P 1 Y 1 
P1P1Q1 



-,n 2 



g(h + b 2 )Z 2 + (3 2 Y 2 _ gb 2 Z 3 - (3 3 Y 3 

, n 3 — 



P2P2Q2 



P3P3Q3 



(15) 



Finally, C,- is determined from the normalization condition 



j{\F^ + \F^)dr = -. 



(16) 



Consequently, we shall gain that G L^R 3 ) at an J t E R and, as a result, 
the solutions of (7) may describe relativistic bound states of quarkonium with 
the energy spectrum (9)— (11). 

Now we can adduce numerical results for constants parametrizing the char- 
monium spectrum which are shown in Table 1. 

One can note that the obtained mass parameter /i is consistent with the 
present-day experimental limits [5] where the current mass of c-quark (2/i ) is 
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Table 2 



Theoretical and experimental charmonium levels. 



C?4- n + n 


3 

Theoret. energy £j = ^2 u>k 
k=i 

(GeV) 


Experim. energy value 
(GeV) 




ei = wi(0, 0, -1) + w 2 (0, 0, -1) + w 3 (0, 0, -1) = 2.979597 


2.979600 


J/V(1S) 


e 2 = Wl (o, 0, -1) + cj 2 (0, 0, 1) + w 3 (0, 0, -1) = 3.096913 


3.096916 


Xco(lP) 


e 3 = wi(0, 0, -1) + w 2 (0, 0, -1) + W3(0, 0, 1) = 3.415186 


3.415190 


Xcl(lP) 


e 4 = wi(0, 0, 1) + w 2 (2, 0, 1) + w 3 (0, 1, -1) = 3.505304 


3.510590 


/lc(lP) 


£ 5 = wi(0, 0, -1) + w 2 (0, 0, 1) + w 3 (0, 0, 1) = 3.532503 


3.526210 


Xc2(lP) 


£ 6 = wi(0, 1, -1) + W2(l, 1, -1) + w 3 (l, 1, -1) = 3.553097 


3.556260 




e 7 = wi(0, 0, 1) + w 2 (l, 0, -1) + w 3 (0, 1, -1) = 3.671608 


3.65400 


Y>(2S) 


e 8 = wi(0, 1, -1) + w 2 (2, 1, 1) + w 3 (l, 1, -1) = 3.674025 


3.685093 


Y>(3770) 


£ 9 = wi(0, 0, 1) + cj 2 (2, 0, -1) + w 3 (0, 0, 1) = 3.775598 


3.770000 


Y>(3836) 


eio = wi(0, 1, -1) + w 2 (0, 0, 1) + w 3 (0, 1, 1) = 3.833640 


3.836000 


X(3872) 


en = wi(0, 1, -1) + kj 2 (0, 1, 1) + w 3 (0, 1, 1) = 3.871672 


3.872000 


Y>(4040) 


£12 = wi(0, 0, 1) + W2(l, 1, 1) + w 3 (0, 1, -1) = 4.042660 


4.040000 


Y>(4160) 


£i 3 = wi(0, 0, -1) + W2(0, 0, -1) + w 3 (0, 1, 1) = 4.153765 


4.159000 


Y>(4415) 


£14 = wi(0, 0, -1) + W2(2, 1, 1) + w 3 (l, 0, -1) = 4.409260 


4.415000 



accepted between 1.1 GeV and 1.4 GeV. As for the gauge coupling constant g 
then its value has been chosen in accordance with many recent considerations 
[6], wherefrom one can conclude that the strong coupling constant a s = g 2 is 
of order 0.22 pa 0.469 2 at the scale of the c-quark current mass. At last, as to 
parameters of solution (3), it is clear that they only shift the origin of 
count for the corresponding energies and we can consider A\ = A 2 = 0. 



With the constants of Table 1 the present-day levels of charmonium spectrum 
were calculated with the help of (9)— (1 1) so Table 2 contains experimental 
values of these levels (from Ref. [5]) and our theoretical ones computed ac- 
cording to the shown combinations (we use the notations of levels from Ref. 
[5])- 
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3 Specification of spectrum: electromagnetic transitions J/ip(lS) 

rj c (lS) +7 and x c0 (lP) — > J ftp (IS) + 7 in dipole approximation 



Now we should specify the obtained above charmonium spectrum. The fact is 
that the relations (9)— (11) permit various parametrizations of the charmonium 
spectrum (see Refs. [2,3]) and therefore it should impose further conditions 
to fix a certain parametrization among several possible ones. In the present 
paper we shall restrict ourselves to elecromagnetic transitions J/ip(lS) — > 
r? c (lS') + 7 and XdoO-P) - ► J/ipO-ty+l since we are interested in the first three 
levels of charmonium. We shall compute widths of the mentioned transitions 
in dipole approximation that will allow us to use the corresponding wave 
functions written out in (13)— (14). Dipole approximation is often employed 
in meson physics, at least, as a primordial width estimate for one or another 
electromagnetic transition (see, e. g. Ref. [7]) so we shall, without going into 
details, remind that in dipole approximation a width T of an electromagnetic 
transition with emission of one photon for some system is given by the relation 
(for more details see, e. g. Ref. [8]) 

r = ^W, (17) 

where dfi - a matrix element of the system dipole moment d = Q(xi+yj+z\c), 
Q is an electromagnetic charge of system, the matrix element is taken between 
a final system stationary state and an initial one and u> is the transition fre- 
quency. Let us introduce the quantities £ = x + iy = r sin^e 1 ^, rj — x — iy — 
r sin^e - ^, z = rcos-#, then it is clear that 

\dfi\ 2 = Q 2 

while in our case Q 2 = (2e/3) 2 with e 2 = 1/137.0359895. We should compute 
the experimental widths T(J/ip(lS) -> rj c (lS) + 7) = r 94 = 1.1557 keV and 
T(xco(lP) -> J/ip( ls ) +7) = r ie = 120.36 keV (data and notations from [5]) 
with the transition frequencies W94 = 0.117316 GeV, cuie = 0.318274 GeV. 



ISM 2 + \vfj\ 2 1 , y ,2 
2 + M 



3.1 Calculation o/r 94 



As is clear from Table 2, the given transition is conditioned by the correspond- 
ing one of the second colour component 0^(0, 0, 1) — > u 2 (0, 0, —1). The initial 
stationary wave function according to (7) is 

* = (19) 
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with the 2D eigenspinor $ 2 °f the euclidean Dirac operator T> 2 on the unit 
sphere § 2 conforming to the eigenvalue A = A 2 = 1 while the final stationary 
wave function is 



■02 — f 1 



r 21^2 



(20) 



with the 2D eigenspinor <3> 2 °f ^2 conforming to the eigenvalue A = A 2 = — 1. 



Further we have F 21 = C 2 P 2 V a 2 e /3 2 r 



1 + 



g(bi+b 2 ) 



with P 2 l = -gih + 62) + P l 2 , Q 



[1 + ^ + 5 2 )] 2 - <? 2 («i + a 2 



?2 

no + u i 2 , Pi 

1 

2 ~~ 



g(fel+&2) 



^-M) 2 + ^ 2 (fei + & 2 ) 2 , 

0^2(0, 0, 1) as well as F 21 = 



C£p£r a *e- I3 * r 



g(bi+b 2 ) 



F f 



1 - 



9(bi+b 2 ) 



with Pi 



- gipi + 6 2 ) + Q2 = Ho + vL Pi = ^l-^if + g^h + hf, a{ = 

\J[—1 + g(Bi + B 2 )} 2 — g 2 (ai + a 2 ) 2 , u> 2 = u; 2 (0, 0, — 1) and we took into ac- 
count that according to (15) Y 2 /Z 2 = —g{b\ + b 2 )/P 2 at n 2 = 0. Also for con- 
stant C 2 we obtain from (16) (with the help of formula [9] / °° x a ~ 1 e~ px dx = 
T(a)p~ a , Rea,p> 0) 



(Cj) s 



1 + 



g{bl+h) )\ptf + (i- g -H^) 2 m 2 



T(2a 2 + 1) _ 1 
P\ J ' \~ P\ J ^ J (2/3*)( 2 <4+D 3' 

(21) 

whereas the analogous relation for C 2 can be obtained from (21) by changing 
indices i — > /. 

At last, as for spinors $ 2 , $ 2 , it should be noted that the eigenvalues of Dirac 
operator T> 2 have the multiplicity 2(/ 2 + 1) for each A 2 = ±(/ 2 + 1), so the 
spinors $ 2 , <£> 2 (l 2 = 0) can be taken in the forms (for more details see Refs. 
[4]) 



S* 1 



c 
2" 



cos 



cos f - i sin f 



c 



cos # — ? sin 
cos I + i sin 2 



|)e^ 2 ,$f = 



C 



cos 2 + i sin | 
2 V — cos I + 2 sin I 



C ( — cos I + i sin I 
cos I + 2 sin I 



(22) 



where the coefficient C = yjl/(2ir). 

Under the circumstances we shall have [(*) signifies complex conjugation] 

R 3 

(23) 

where ^(^2)^2 * s the sum over a ^ possible combinations of final and initial 
spinors. It is then not complicated to see that the only nonzero contribution 
comes from the combination ($ 2 1 )''"$ 2 2 = e~ tLp {iC 2 sini?)/2 while the rest of 
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combinations give zero when integrating over tp in (23). Fulfiling the remaining 
integration in (23) over i? and r yields 



2 



i— C\C 2 



T(a f 2 +a i 2 + 2) 
(Pi + P l 2 ) a * +a * +2 



x 



i + 



gih + & 2 ) 

Pi 



i + 



^(&i+& 2 ) 



+ Q2Q2 



gih + b 2 ) 

Pi 



g(h + b 2 ) 



p\ 

(24) 

The similar considerations show that rjfi ^ only for the spinor combina- 
tion ($^ 2 )" i "$2 1 = — e iV (iC 2 sin$)/2 which entails £/j = —rjfi while two spinor 
combinations giving the contributions to z/i are (& 2 2 y&f — ~ (^^V^i = 
cos??/ (87r) which gives rise to z/i = 0. As a result, we obtain according to (18) 



— 



4 ^4 ^2 



(25) 



with £ fi of (24). 



3.2 Calculation o/T 16 



In accordance with Table 2 the given transition is formed by two ones: the 
photon absorption with frequency oji = cj 2 (0,0, 1) — u; 2 (0,0, —1) = 0.117316 
GeV by the second colour component and the photon emission with frequency 
10 11 = 003(0, 0, 1) — ^3(0, 0, — 1) = 0.435591 GeV by the third colour component 
so the resulting frequency ujiq = ion — uj = 0.318275 GeV. Accordingly we 
can write T 16 = T/ + F n . Now the considerations analogous to the above ones 
show that again Cfi — ~Vfii z fi — m both transitions and Ti is obtained 
from (24)-(25) by replacing % —>/,/—> i, ^94 — > while Y n will be the 
result of changing indices 2^3 and {pi + b 2 ) — > — b 2 , (a\ + a 2 ) — > — a 2 , 
(Si + B 2 ) -> -S 2 , CJ94 -> wj/ in (24)-(25). 



3.3 Numerical results 



Table 3 gives the results of evaluation for r 94 and T 16 according to the ob- 
tained relations and the corresponding experimental ones of Ref. [5]. Table 
3 supplies us with an additional justification for the choice of parameters of 
the SU(3)-confining gluon field adduced in Table 1 and allows us to conclude 
that dipole approximation is not enough for the second transition of the ones 
under consideration. The question now is what gluon concentrations are in 
the mentioned SU(3)-confining gluon field. 
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Table 3 

Widths of electromagnetic transitions in charmonium. 



r 


Theoret. value 


Experim. value 




(keV) 


(keV) 




1.02977 


1.15570 


Tie 


6.06961 


120.360 



4 Estimates of gluon concentrations 



To obtain necessary estimates we shall use an analogy with classical elec- 
trodynamics where is well known (see e. g. [10]) that the notion of classical 
electromagnetic field (a photon condensate) generated by a charged particle 
is applicable only at distances >> the Compton wavelength A c = 1/m for 
the given particle. If denoting Xb the de Broglie wavelength of the particle 
then Xb = 1/p with the relativistic impulse p = mv/yl — v 2 while v is the 
velocity of the particle (as a result, A c = Xb at v = l/\/2) so one can rewrite 
Xb = X c \/1 — v 2 /v and it is clear that Xb — > when v — > 1 (ultrarelativistic 
case), i. e., the particle becomes more and more point-like one. Accordingly, 
one can conclude that in the latter case the notion of classical electromagnetic 
field generated by a charged ultrarelativistic particle is applicable at distances 
>> Xb- Under the circumstances, if the ultrarelativistic charged particle ac- 
complishes its motion within the region with characteristic size of order r 
then in the given region the electromagnetic field generated by the particle 
may be considered as classical one at r >> Xb- For example, in the case of 
positronium we have tq ~ 2ao, where ao ~ 5.29 • 10~ n m is the Bohr radius, 
so r >> A e ~ 3.86 • 10~ 13 m, the electron Compton wavelength, i. e., the 
electric Coulomb interaction between electron and positron in positronium can 
be considered as classical electromagnetic field. 

Passing on to QCD, gluons and quarkonia, it should be noted that quarks in 
quarkonia accomplish a finite motion within a region of order 1 fm = 10~ 15 m. 
Then, as is seen from the radial parts of the wave functions (13)-(14), the 
quantity l/j3j permits to be considered as a characteristic size of the j-th 
colour component of the given quarkonium state and, consequently, we can 
take the magnitude 

r °=lU < 26) 

for a characteristic size of the whole quarkonium state and, in line with the 
above, we should consider the confining SU(3)-gluonic Yang-Mills field of (3) 
or (6) to be classical one when r >> X B , the de Broglie wavelength of the 
corresponding quarks forming quarkonium. 

On the other hand, a classical electromagnetic field (photon condensate) con- 
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forms to the large photon concentrations for every frequency presented in the 
field [8] . Then in QCD we should require the large gluon concentrations in the 
given classical gluonic field (gluon condensate). To estimate the given concen- 
trations we can employ T 00 -component of the energy-momentum tensor of (5) 
and, taking the quantity u = T, the whole decay width of the quarkonium 
state, for the characteristic frequency we obtain the sought characteristic con- 
centration n in the form 

Too / _x 

It is not complicated to obtain the curvature matrix (field strentgh) corre- 
sponding to the solution (3) or (6) 

F = F^Xadx" A dx u = -d r {A a t \ a )dt Adr + d r (A a ^X a )dr A dip , (28) 

which entails the only nonzero components 

3 _ 2a 1 + a 2 8 _ aW3_ 3 _ h 8 _ jW3_ 
r tr - ' tr ~ 2r 2 ' 9 ~ 2 9 ~ 2 

and, in its turn, 



3 (aj + a x a 2 + ai b[ + b x b 2 + b 2 2 \ 



16n \ r 4 3r 2 sin 2 

so, further putting sin 2 ?? = 1/3 for simplicity, we can rewrite (30) in the form 

T 00 = T c o oul + Tqq 11 (31) 

conforming to the contributions from the Coulomb and linear parts of the 
solution (3) or (6). The latter gives the corresponding split of n from (27) 

n = n coul + n ha . (32) 

For comparison we shall also estimate the photon concentration in the ground 
state of the positronium. As is known historically [11], the analogy between 
the latter system and quarkonia played the important part when building 
the quarkonia models. For positronium we have the electromagnetic Coulomb 
interaction A = A t dt = (e/r)dt which entails F = F tr dt A dr = (e/r 2 )dt A dr 
and 

Too = T U = l(F lr f = ^ (33) 
with a em = e 2 = 1/137.0359895. 



11 



5 Numerical results 

When computing for the first three states of charmonium we used their present- 
day whole decay widths T = 17.3 MeV, 91.0 keV, 10.2 MeV respectively [5], 
while the calculation r of (26) gives r = r 1 ,r 2 ,r 3 (see Table 4). 

In the positronium case we employed the widths T = l/r (parapositronium) 
and Fi = 1/ti (orthopositronium), respectively, with the life times t = 1.252- 
10~ 10 s, Ti = 1.377 • 10~ 7 s [11] while r = 2a with the Bohr radius a = 
0.529177249- 10 5 fm [5]. 

Tables 4, 5 contain the numerical results for both the cases, where we used 
1 A = 10~ 10 m = 10 5 fm. Finally, when calculating we applied the relation 
1 GeV- 1 « 0.21030893 fm . 



6 Concluding remarks 

It is clear that for charmonium the natural unit of the gluon concentration 
is fm -3 while for positronium the photon concentration should be measured 
in A . Then, as is seen from Tables 4, 5, qualitative behaviour of both the 
concentrations is similar. At the characteristic scales of each system the con- 
centrations are large and the corresponding fields (electric and magnetic colour 
ones or electric Coulomb one) can be considered as classical ones. For char- 
monium the part n cou i of gluon concentration n connected with the Coulomb 
colour electric field is decreasing faster than nn n , the part of n related to the 
linear colour magnetic field, and at large distances n\ in becomes dominant. 
Under the circumstances, as has been said in Section 4, we can estimate the 
quark velocities in the charmonium states under discussion from the condition 



with the c-quark Compton wavelength \ q = l/(2/i ) 0.168247 fm so, taking 
the de Broglie wavelength A# = O.lr^ (k = 1,2,3) with from Table 4, 
we obtain v x « 0.999718, v 2 « 0.999721, v 3 « 0.999715 i. e., the quarks in 
charmonium should be considered as the ultrarelativistic point-like particles. 
This additionally confirms the conclusion of Refs. [2,3] that the relativistic 
effects are extremely important for the confinement mechanism. As a result, 
the confinement scenario described early in Section 1 may really occur. 

The results of the given paper can be extended to both other states of char- 
monium and the bottomonium states. We hope to return to the problem else- 
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Table 4 

Gluon concentrations in the first three states of charmonium. 



»7 C (1S): ri = 0.0399766 fm 



r n cou i nii n n 

fm (far 3 ) (fm" 3 ) (far 3 ) 



O.ln 0.550649 • 10 13 0.727630 • 10 iU 0.551377 • 10 13 

n 0.550649 • 10 9 0.727630 • 10 8 0.623412 • 10 9 

10ri 0.550649 • 10 5 0.727630 • 10 6 0.782695 • 10 6 

1.0 0.140637 • 10 4 0.116285 • 10 6 0.117691 • 10 6 

a 0.179347 • 10~ 15 0.415260 • 10~ 4 0.415260 • 10~ 4 



J/^(15): r 2 = 0.0397797 fm 



r n cou i nii n n 

fm (fm" 3 ) (fm -3 ) (far 3 ) 



0.1r 2 0.106772 • 10 lb 0.139702 • 10 i3 0.106912 • 10 lb 

r 2 0.106772 • 10 12 0.139702 • 10 11 0.120742 • 10 12 

10r 2 0.106772 -10 s 0.139702 • 10 9 0.150380 • 10 9 

1.0 0.267364 • 10 6 0.221069 • 10 8 0.223742 • 10 8 

a 0.340956 • 10~ 13 0.789450 • 10~ 2 0.789450 • 10~ 2 



Xco(l-P): r 3 = 0.0401851 fm 



r n cou \ n un n 

fm (far 3 ) (fm -3 ) (fm" 3 ) 



0.1r 3 0.914712 • 10 13 0.122134 • 10 11 0.915933 • 10 13 

r 3 0.914712 • 10 9 0.122134 • 10 9 0.103685 • 10 10 

10r 3 0.914712 • 10 5 0.122134 • 10 7 0.131282 • 10 7 

1.0 0.238531 • 10 4 0.197228 • 10 6 0.199613 • 10 6 

a 0.304186 • 10~ 15 X ?).704313 • 10~ 4 0.704313 • 10~ 4 



Table 5 

Photon concentrations in the ground state of positronium. 



r = 2a = 2 • 0.529177249 • 10 5 fm 



Par ap ositronium 



Orthopositronium 



r 



para 




fm 



(fm"' 




(fm- 



0.01r 0.888025 • 10 (0.888025 • 10 16 ) 0.976685 • 10 4 (0.976685 • 10 19 ) 



r 0.888025 • 10~ 7 (0.888025 • 10 8 ) 0.976685 • 10~ 4 (0.976685 • 10 11 ) 
2r 0.555015 • 10~ 8 (0.555015 • 10 7 ) 0.610428 • 10~ 5 (0.610428 • 10 10 ) 



where. 
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